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A vortex crossing a thin-film superconducting strip from one edge to the other, perpendicular to the bias 
cun'ent, is the dominant mechanism of dissipation for films of thickness d on the order of the coherence length 
^ and of width w much narrower than the Pearl length A 2> w ^ C- ^.t high bias currents, I* < I < Ic, 
the heat released by the crossing of a single vortex suffices to create a belt-like normal-state region across the 
strip, resulting in a detectable voltage pulse. Here Ic is the critical current at which the energy barrier vanishes 
for a single vortex crossing. The belt forms along the vortex path and causes a transition of the entire strip into 
the normal state. We estimate J* to be roughly /c/3. Further, we argue that such "hot" vortex crossings are 
the origin of dark counts in photon detectors, which operate in the regime of metastable superconductivity at 
currents between /* and Ic- We estimate the rate of vortex crossings and compare it with recent experimental 
data for dark counts. For currents below J*, i.e., in the stable superconducting but resistive regime, we estimate 
the amplitude and duration of voltage pulses induced by a single vortex crossing. 



I. INTRODUCTION 

Dissipation in superconducting wires thinner than the 
coherence length ^ have been thoroughly studied both 
theoreticall}^ and experimentally;^ In these one-dimensional 
(ID) superconductors the dissipation arises due to 27r-phase 
slips occurring in segments of length ^ of a wire that becomes 
temporarily normal. Langer and Ambegaokar— treated the 
problem of dissipation in ID wires with ring geometry within 
the theory of nucleation rates of current-reducing fluctuations 
in a superconductor. The transition between states with dif- 
ferent currents in a ring occurs via the nonstationary state de- 
scribed by the saddle point solution of the Ginzburg-Landau 
(GL) functional. Langer and Ambegaokar' found such a so- 
lution and the corresponding free energy difference or bar- 
rier, U, between the original metastable state with current and 
the saddle point state (see also Ref. 0)- Later McCumber and 
Halperin derived the attempt frequency ft in the phase-slip 
rate, R = ilcxp{—U /T), using time-dependent GL theory.— 

The problem of dissipation in superconducting thin-film 
strips with the thickness d much smaller than the London pen- 
etration depth A, and of width w much smaller than the Pearl 
length, A = 2X^/d ^ w, has been extensively discussed 
in the context of a possible Berezinsky-Kosterlitz-Thouless 
(BKT) transition in superconducting filmsi^— The interest in 
current-carrying thin-film strips has been revived recently in 
search for quantum tunneling of vortices.r'— their dynamic 
behavior^ia and the observation of so-called "dark counts" in 
superconductor-based photon detectors. '^- '"^ The detector con- 
sists of a long and thin superconducting strip carrying cuiTents 
slightly below the critical value. Typically, in NbN photon de- 
tectors w is of the order of 100 nm or more and c? 4 — 6 nm, 
while the zero-temperature coherence length ^(0) ss 4 nm. 
The low-temperature London penetration depth A « 350 nm 
so that the Pearl lengtbi^ A w 40 /im > w. 

When a photon interacts with the strip it induces a hot spot 
in the film that drives a belt-like region across the strip in the 
normal state. Consequently, a voltage pulse caused by the 



current redistribution between the superconducting strip and a 
parallel shunt resistor is detected. After the normal belt of the 
strip cools down, the strip returns to superconducting state. 
Thus, single photons can be detected and counted by measur- 
ing voltage pulses. However, similar pulses are recorded even 
without photons (dark counts). These voltage pulses have 
peak amplitudes similar to photon-induced pulses."' There- 
fore, one can conclude that dark counts are also caused by 
nucleation of normal belts across the strip. In both cases and 
in the absence of a shunt, the entire strip undergoes transition 
into the normal state due to heat released by the bias current 
in the normal belt region. 

In fact, the observation of dark counts means that the su- 
perconducting strip, at bias currents slightly below the critical 
current, is in a metastable state. Photons or fluctuations trigger 
the transition from this state to the normal state. Thus, the cen- 
tral question is what kind of fluctuations trigger the transition 
in the case of dark counts. The origin of dark counts is still 
debated (see Refs. [Tsl and [13). The problem of dark counts 
is related to the basic question of dissipation in thin films and 
wires and is of technological relevance because fluctuations 
resulting in the formation of normal belt across the strip limit 
the ability of superconducting circuits to carry supercurrents, 
in general, and the accuracy of photon detectors, in particu- 
lar. In the literature, dark counts are treated either within the 
formal framework of ID phase slips in thin wires or within 
the picture of vortex-antivortex unbinding near the BKT tran- 
sition (see Refs. ITsI andlTil). Vortices crossing the strip were 
employed to explain dc current-voltage characteristics of thin- 
fikn stripsjiiiii 

In this paper we discuss three types of possible fluctuations 
in superconducting strips which result in dissipation. Each 
one causes transition to the normal state from the metastable 
superconducting state when currents are close to the critical 
value Ic- 

(a) Spontaneous nucleation of a normal-state belt across the 
strip with 27r-phase slip as in thin wires. 

(b) Spontaneous nucleation of a single vortex near the edge of 
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the strip and its motion across to the opposite edge accompa- 
nied by a vohage pulse. 

(c) Spontaneous nucleation of vortex-antivortex pairs and 
their unbinding as they move across the strip to opposite edges 
due to the Lorentz force, as well as the opposite process of nu- 
cleation of vortices and antivortices at the opposite edges and 
their annihilation in the strip middle. 

The energy barrier for the nucleation of a temporary normal 
phase-slip belt is too high to be of importance because the belt 
volume > d^w is large. We will show that such a barrier re- 
mains large at any current in the superconducting state. Con- 
sequently, belt-like 27r-phase slips appear with extremely low 
probability. On the other hand, as proposed in Refs.fTTlandflTl 
thermally induced vortex crossings in current-carrying strips 
result in 27r-phase changes along the strip just as in the ID 
scenario and hence cause dissipation. For the case of quan- 
tum tunneling this mechanism of dissipation was discussed 
in Refs. [st-UH The free energy barrier for vortex crossing is 
much lower than for belt-like 27r-phase slips, since the vor- 
tex core volume is d£^^ <C d^w. The energy cost of creating 
a vortex and moving it over the barrier is uj/^ times smaller 
than for creating a belt-like phase slip. An important point is 
that such a barrier for vortex crossing vanishes as the current 
approaches Ic, whereas the barrier for the belt-like phase slip 
remains nonzero at any current. As to the vortex-antivortex 
process of the point (c), we show in the following that the 
corresponding barrier is twice as high as for the single vortex 
process. 

We evaluate the amplitude of a voltage pulse and its dura- 
tion assuming that the belt-like area around the vortex path 
remains superconducting. We call this process a "cold" pulse. 
This is not always the case, because vortex motion excites 
quasiparticles along the vortex path and their energies depend- 
ing on the bias current may suffice for creation of a normal- 
state belt across the strip. This will result in redistribution 
of current from the superconducting strip to the shunt with 
the accompanied voltage pulse much bigger than for "cold" 
pulses. Such a "hot" pulse will be similar to the one induced 
by photons. In the following we will estimate at what min- 
imum bias current /* a single vortex crossing can trigger a 
"hot" voltage pulse and a corresponding dark count. 

Thus, we argue that dissipation and corresponding voltage 
pulses in strips are caused predominantly by vortex crossings. 
At high bias currents such crossings release energy sufficient 
for the formation of a normal belt along the vortex trajectory, 
see Fig. \lla). Such a belt triggers the transition of the whole 
strip into the normal state in the absence of a shunt resistor, as 
well as the redistribution of the bias current into the shunt in 
the case of photon detectors. Note that a similar process hap- 
pens when a photon creates a normal "hot" spot on the strip. 
When this spot is sufficiently large, it destroys the supercon- 
ducting path for the transport current and the current redistri- 
bution leads to a voltage pulse, the photon count. If the hot 
spot does not disrupt completely the superconducting path, it 
will nevertheless lead to a decrease of the energy barrier for 
subsequent vortex crossings. At high bias currents, a "hot" 
vortex crossing can happen directly, see Fig. [Ha), or through 
a hot spot area created by photon and forming a normal belt. 
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FIG. 1 . (Color online) Sketch of a segment of the strip in the presence 
of a bias cun'ent /. Panel (a): a single vortex (blue circle) causes a 
"hot" crossing (pink belt). The width of the belt £ is of the order of 
superconducting coherence length. Panel (b): a single photon creates 
a hotspot (red disc) and induces a subsequent "hot" vortex crossing 
(pink belt). Both processes result in detectable voltage pulses in a 
superconducting nanowire single photon detector (SNSPD). 



which will result in signal detection, see Fig.fTJb). 

The layout of this paper is as follows: In Sec. II we discuss 
three energy barrier scenarios for vortex crossings. In Sec. Ill 
we derive dc current-voltage characteristics and evaluate the 
magnitude of induced voltage pulses. The concept of "cold" 
and "hot" vortex crossings is introduced in Sec. IV. In Sec. V 
we compare our results with data for dark count rates in NbN 
filmsji^ We summarize our results in Sec. VI. 



II. ENERGY BARRIERS AND VORTEX CROSSINGS 

In this section we derive energy barriers for three dissipa- 
tive processes mentioned within the GL theory. Consider a 
thin-film strip of width w A and of length L ^ w. We 
choose the coordinates so that < x < w and —L/2 < 
y < L/2. Since we are interested in bias currents which may 
approach depairing values, the suppression of the supercon- 
ducting order parameter must be taken into account. We use 
the standard GL functional with respect to the order parame- 
ter 4'(r) (normalized to its zero-field value in the absence of 
current) and the vector potential A: 
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Here $o is the flux quantum, r ~ (x, y) is a point on the 
film, V is the 2D gradient, and He = $o /2\/27tX^ is the ther- 
modynamic critical field. The order parameter in the presence 
of a uniform bias current / in zero applied magnetic fields and 
with no vortices present, can be found by minimizing the GL 
functional and disregarding the current self-field, as is done, 
e.g., in Ref.0- As discussed in the next section, this is an accu- 



3 



rate approximation for w A. Thus we obtain the solution: 

= (1 - K2)l/2g-z,cy/C+^<po^ (2) 



where (p^ is an arbitrary constant phase. The parameter k is 
proportional to the phase gradient and describes the order pa- 
rameter suppression due to bias current. As a function of k 
the bias current in the superconducting state is limited to the 
depairing current /max = /o(4/37r-\/3)('u;/^), corresponding 
to Kmax = l/'v/S, as for the case of the ID wire. 



A. Phase slip in the normal belt 

When dealing with the situation of fixed uniform current / 
instead of vector potential A it is more suitable to work with 
the Gibbs free energy functional rather than the free energy 
functional, Eq. ([T]). We perform the usual Legendre transform 
(see Ref. [H) to obtain the corresponding free energy density: 
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The equilibrium Gibbs free energy density for a given current 
is obtained by minimization with respect to *. It jumps at the 
maximum current /,nax from //(/max) = ~(2/9)(//^/87r) to 
zero at /max as expected for a first order transition. Hence, the 
free energy barrier U for creation of a belt-like normal-state 
area with volume V = iwd (£ is the width of the belt along 
the y-axis) decreases from {H^/8n)V to {2/9){H^ /8tt)V as 
the bias current increases from to /max- The barrier never 
vanishes in this interval ("overheating" with respect to bias 
current is absent). Note that for w ^ £, and £ > £, the barrier 
remains very high in comparison with the temperature at all 
bias currents I < Ic resulting in low probability for phase 
slips, except for temperatures close to Tc, where the barrier 
vanishes as (1 - T/Tc)^. 



B. Single vortex crossing 

A vortex crossing from one strip edge to the opposite one 
induces a phase slip without creating a normal region across 
the strip width. We will treat the vortex as a particle moving 
in the energy potential formed by the superconducting cur- 
rents around vortex center inside the strip and by the Lorentz 
force induced by the bias current. We will derive the energy 
potential and find the vortex crossings rate (phase slips and 
corresponding voltage pulses) in the framework of Langevin 
equation for viscous vortex motion and invoke the known so- 
lution of the corresponding Fokker-Planck equation. 

In the presence of a vortex, the order parameter in the 
current-carrying strip, disregarding its suppression in the vor- 
tex core, reads: 



(5) 
(6) 



In this approximation, the vortex affects mainly the phase 
if{r,r„) of the order parameter To describe voltage pulses 
we need to know how the phase changes when the vortex 
moves across the strip. For simplicity we consider a vortex 
at {xvjUv) = 0. As we ignore the change of the order pa- 
rameter amplitude in the vortex core, the current distribution 
is governed by the London equation (integrated over the film 
thickness): 



+ 27r(A/c) curiae/ = $o ^(r- - r„), 



(7) 



where g is the sheet current density. 

For narrow strips, w <C A, the field is approximately 
hz ^ g/c, whereas the term with derivatives is of the order 
Ag I cw. Hence, in this limit, supercurrents can be found by 
neglecting and the corresponding vector potential of the 
order w/A,i2i2S. Introducing the scalar stream function G(r) 
such that 

g = curl(Gi) (8) 
we reduce the problem to solving the Poisson equation: 



V^G = -(c$o/27rA)5(r - r„) . 



(9) 



Since the boundary condition at the strip edges requires 
vanishing normal components of the current, we have G = 
at X = 0, w. Therefore, the problem is equivalent to one in 
2D electrostatics: a linear charge at r„ between two parallel 
grounded plates atx = 0,w with the known solution)^ 

G(r) - In '^Qshy - cos{X + Xy) 
^ n '\oshy-cos(X- A:^,) ' 

where capitals stand for coordinates in units of w/tt, i.e., x = 
X w/tt, y ^ Y wjix. 

The energy of a vortex at a; = .Tu and y = is: 



$0 

2c 



G(a;,,0), 



(10) 



with the standard cutoff ^ at the vortex core.— In the presence 
of a uniform bias current the energy barrier reads: 



U(X,) 
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sin X, 
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(11) 
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where eq is the characteristic energy of a vortex in thin 
films. The vortex energy U{Xj,) is maximum at Xg ~ 
tan~ ^ (h^Iq 1 1) and the energy barrier is given by 
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■ tan 
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(13) 

This barrier decreases with increasing current and turns zero 
at a critical value on the order of the depairing GL current: 



/. 



■ wd. 
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here e = 2.718. One can see that the critical current Ic is 
slightly smaller than /max discussed above. 

Since the vortex mass is negligibly small, we use the equa- 
tion of purely diffusive motion (only includes first order time 
derivative) for describing the vortex propagation between x = 
and X ^ w: 



dX,, 



7- 



dt 



where 7 = vo^rj/iT^ and 



dUiX,) 
dX^ 



^0 



(15) 



(16) 



is the Bardeen-Stephen drag coefficient for fibn with = 
Pn/d being the film's sheet resistance slightly above Tc- F{t) 
is the Langevin random force obeying statistical averages 
{F{t)) ^ and {F{t)F{t')) = 2jTSit - t'). 

The vortex motion described by Eq. ( fTST i, takes place in the 
interval a < x < w — a, where a is of the order of ^ (the en- 
ergy of the system cannot be described by the potential (fTTl i in 
the intervals w — a < x < w and < a; < a). The most cru- 
cial interval for vortex motion is near the point Xs = Xsw/tt, 
where the vortex should overcome the potential barrier. Thus 
Xs should be inside the interval (a, w ~ a), i.e., the conditions 
^ <^ w and / < (e/2)/c should be fulfilled to consider the 
motion of vortex in the interval < x < u;. To compute the 
average velocity in the interval < x < w, we consider the 
diffusion problem of a single particle that propagates in the in- 
terval —00 < a: < 00 under the effect of the periodic potential 
ev{x) ~ ey{x + w) and the Lorentz force. The average ve- 
locity is obtained from the known stationary solution for this 
periodic model (see Ref. I22I) . This approach was previously 
used by Gurevich and Viinokur— . 

The corresponding Fokker-Planck equation (Smoluchowski 
equation) for the probability current in the case of the periodic 
potential has a stationary solution with the statistical average 
vortex velocity v given hy^ 

^jj^ ^ (17) 

Z±{x)= f due^^^"'^''^, P^l-e-'^P (18) 

where v = X and p = i/I/p^Iq. Except for temperatures 
close to Tc the parameter = ^^eq/T ^ 1 . At large v the 
function exp^{x)/T] has a sharp maximum between and 
w, while the function cxp[—U{x)/T] has two sharp maxima 
at the edges of this interval. Since the integral Z+ (tt) has the 
analytic solutions^ 



dxe P^sin"; 



7rexp(-7rp/2)r(i^ + 1) 
2''\T{l + iy/2 + ip/2)\^' 



(19) 



where r(a:) is the Gamma-function and ly > —1, we obtain 
the asymptotic solution for » 1: 



2w 



P 



-p tan 



(20) 



Evaluating Z_ (tt) we note that the main contribution comes 
from the regions near the edges, where we approximate 
sin(j:) ~ sm{n — x) ^ x and replace the low integration 
limit by TT^/ w and the upper one by tt — tt^/ w. We obtain the 
asymptotic limit 

(21) 
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In the integral /J" da;e~"(^)/^Z+(.T), the function Z+{x) 
reaches maximum at x = n and is small at low x. Hence, 
the main contribution to this integral comes from the region 
near x ^ tt: 

'1 



1 



Z+{n). 
(22) 

It then follows that the dependence of the average vortex ve- 
locity w on / at large p and v is given by 

Y{z) = (l + z2)('^+i)/2cxp[zyztaii-i(l/z)]. (24) 
Note the strong power-law dependence of v on the strip width 



For large currents, I ^ Iq, this expression reduces to 



7ti 



T 



(25) 



2 J (<) (4, 
with Ic given by Eq. (fT4] i. Note that the average velocity 
changes drastically near the critical current Ic, where the en- 
ergy barrier vanishes. Such defined critical current is about 
16% smaller than the standard depairing current /max defined 
for 1 D wires (vanishing energy barrier for phase slips in wires, 
see Ref.i),24 

In the case of multiple simultaneous vortex crossings hap- 
pening in different parts of the strip, we must account for their 
interactions. The interaction of vortices situated at (Xi,0) 
and {X2, Y) has been evaluated in Ref.O 



eint = eo In 



coshF - cos(Xi +X2) 
coshy — cos(Xi — X2) 



(26) 



If vortices are separated hy y > w along the strip, the in- 
teraction is exponentially weak and their crossings are uncor- 
rected. Accounting for both vortex and antivortex crossings 
(which are equivalent by symmetry), we estimate the rate for 
multiple vortex crossings at / < /^ as /? ss {2L/'kw)v. 
Finally, we obtain the asymptotic estimate for the rate: 



R: 



TTiy 



1/2 



y 



I 



(27) 



In obtaining this result we disregarded vortices crossing in 
the direction opposite to the Lorentz force, the correspond- 
ing probability for such processes is oc e^'^P <C 1. We note 
that Gurevich and Vinokur took L/^ as the number of statis- 
tically independent vortex crossingS/i2, It differs by a factor 
^/w <^ 1 from our estimated number L/w of independent 
crossings. Therefore, Ref . (Tt] overestimates the rate. 
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C. Vortex-antivortex pair scenario 



The energy of a vortex-antivortex pair (vortex-antivortex in- 
teraction included) was derived in Ref. [T9| and is 



In 



4PF2 coshy-cos(Xi -X2) 

sm Ai sinA2 



2f2 



coshy - cos(Xi + X2) 



(28) 

This energy increases with increasing separation Y, so that 
one expects the lowest barriers for Y ^ 0: 



= In 



sin Xi sin X2 



sn-?[{Xi-X2)/2] 



'■[{Xi+X2)/2] 



(29) 

One can show that if a pair is formed at Xq and the pair mem- 
bers are pushed apart a distance 26, the lowest energy increase 
(for a given b) corresponds to the initial position Xo = 7r/2 
in the middle of the strip. The energy barrier for such a 
pair, in the presence of bias current /, is obtained by setting 
Xi^2 ~ 7r/2 =F h and adding the Lorentz force contribution: 



In 



w sin(26) 



lb 

Wo 



(30) 



This energy is maximum if 2b = tan ^ (2/1^/0 / /) so that the 
energy barrier for vortex-antivortex pairs is given by 



Hp 



-In 



1 



■ tan 



(31) 

For I ^ Iq this barrier is twice as large than that for a sin- 
gle vortex crossing, Eq. (fT3] l, and the ratio of these barriers 
increases for smaller currents. Note also that the core contri- 
bution to the pair energy (neglected here) is at least twice that 
for a single vortex. 

Based on our estimates for the three different fluctua- 
tion scenarios presented here, we conclude that single vortex 
crossings are the main source for dark counts. 



III. VOLTAGE INDUCED BY VORTEX CROSSING 

Let us now find how the phase of the order parameter varies 
when a vortex crosses the strip. The current is expressed ei- 
ther in terms of the gauge invariant phase ip or via the stream 
function G: g = -{c^o/4:Tr^A)S/Lp = curl[G'£]. Written 
in components, this gives the Cauchy-Riemann relations for 
functions [Att'^ Aq / c^Qfi'^]G{r) and (p{r). Hence they are real 
and imaginary parts of an analytic function of complex argu- 



ment z = X - 



ly 



.21 



g{z) = hi 



sin[(X„ + Z)/2] 
s\tl[{X^ - Z)/2] 



(32) 



(recall: the capitals are coordinates in units of w/tt, so that 
< X < TT, etc.). We then obtain 



^(r,r„) =Im[g(Z) 
— tan 



(33) 



s\nXy sinh(y — Y^) 



cos X — cosh(F — Ytj) cos X^ 



Note that the characteristic length of variations for Lp in both 
X and y directions is w. For long strips of interest, L ^ w, 
and for distances |y — Yuj ^ 1, we have at the strip ends 
Lp{±L/2) = ^Xy. Hence, when the vortex moves from the 
strip edge at Xy = to the opposite edge at Xy = n and 
\L/2 — Yy\ 1, the phase difference at the ends of the strip 
changes by ip{L/2) — ip{—L/2) = 2Xy = 27r, i.e., a vortex 
crossing results in a global phase sUp of 27r. 



A. DC voltage 

The motion of vortices causes the phase difference at the 
strip ends to vary in time. Using the Josephson relation for 
the phase, we obtain the induced voltage due to a single vortex 
crossing 



(34) 



where the vortex velocity is v{t) = dxy/dt — {w/TT)dXy/dt 
and we used Lp{L/2) — ip{—L/2) = 2Xy. A quasistatic ap- 
proach employed here is justified as long as the characteristic 
crossing time At ^ w/v is large compared to L/c. Note that 
for each crossing, i.e., for each voltage pulse between time t 
and t + At the relation 



t+At 



dt' V{t') = ^ 



(35) 



is satisfied as in the case of voltage pulses due to phase slips 
in ID wires.* Thus we obtain the average (dc) voltage 



$0 



R. 



(36) 



This relation also follows directly from comparing the dissi- 
pated power Vdc^ with the work per unit time done by the 
Lorentz force, {^qI / cw)wR. It is worth to remember that we 
have derived the crossing rate assuming an isothermal strip . In 
continuous measurements of current-voltage characteristics at 
currents of the order of the critical one, the strip temperature is 
certainly higher than that of the bath. In principle, this heating 
may be reduced using short bias current pulses. 



B. Voltage pulses 

In this section we consider the time evolution of the voltage 
pulse V{t) induced by single vortex crossing. Here we use the 
equation of vortex motion, Eq. (flST l. for X > Xg and neglect 
random forces (thermal noise). Therefore the velocity is 



— [—- H^coiX 

T]W \Io 



This can be written in the form 



X = /3(cotX^-cotX), /3 
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(37) 



(38) 



which is vaUd for X > Xg- Itis worth noting that for currents 
of the order of the saddle point is very close to the strip 
edge. 
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Integration of Eq. 



results in an implicit solution for: X{t): 

X (t) cos X, + sin Xs In sin [X (t) - X J = ^ "J"^ . 

sm Ac, 

(40) 

We choose the constant to so that t = corresponds to 
the vortex exit at X = tt. Note that any instant for which 
< X{t) < Xs is beyond this approximation, because in 
this early time interval the process is described by thermal 
activation rather than by the equation of motion (|37] | with ran- 
dom force omitted. The instant for which X{t) = Xs is also 
inappropriate as an initial moment, because at this point the 
velocity vanishes, X ^ 0. Thus Eq. ( l40b can be written as 



[X{t) - tt] cos Xs + sin X, In 



sin [Xit)-Xs] 
s'mX^ 



f3t 
sin Xs 
(41) 



Clearly, X{0) = tt and X{t — oo) = Xs- Hence, formally, 
the motion from the saddle point Xg to the edge takes infi- 
nite time because the velocity goes to zero as X ^ Xs- In 
reality, the dynamic viscous vortex motion starts at some dis- 
tance from the saddle point where the vortex is kicked by ran- 
dom force (an activation driven process) and the total "time- 
of-flight" is finite. To see this, consider the situation of large 
currents for which Xs is given by Eq. ( [39] l and 



X{t) 



Xs\n 



sin [X{t)~Xs] 



f3t 

'X's' 



(42) 



smXs 

Denote as 5X a small distance from the saddle at Xs and eval- 
uate the time tq of motion from Xg + 5X to the edge X = tt: 



Xs 



Xs + 6X 



Xs In ■ 



sx_ 
XT' 



Since both SX and Xs are small, all terms on the right hand 
side, except for tt, are negligible and we obtain 



To 



ttXs 



crjw 



(43) 



so that the time-of-flight tq does not depend on a particular 
choice of 5X- In fact, this estimate coincides with the time 
it takes a vortex to cross the strip being pushed solely by the 
Lorentz force. 

Solving numerically Eq. (l4Tl i for X{t) and substituting the 
result in Eq. ( l38T l we obtain v{t)- The result is shown in Fig.|2] 
For convenience, we use Xs/P as the unit of time. The dimen- 
sionless time t = Pt/ Xs varies between — tt < t < 0. 

The divergence at the edge x = w must be cut off at 
distances of the order of ^ from the edge. We obtain from 
Eq. (1381 ) an estimate for the maximum velocity at the exit. 
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FIG. 2. (Color online) The dimensionless vortex velocity dX/dr 
versus r for parameters Xs = 0.1, 0.3, and 1.0. Note that dX/dr = 
v/v, where v is the velocity in common units and v is the average 
velocity, which is identical to the one solely due to the Lorentz force. 



where the critical current is given by Eq. (Il4i . 

For large currents, Xs ^ 1, we solve Eq. (l42l l perturba- 
tively: X ^ Xi+5X with Xi = tt + pt/Xg and 5X X-^: 



X 



T - Xg In 



sin(X, 



sinXs 



Thus the velocity for Xg <C 1 is 

dX 



^1+XsCOtiXs-T). 



(45) 



(46) 



the unity corresponds to a constant velocity due to the Lorentz 
force, whereas the second term is caused by the vortex poten- 
tial. 

The velocity v{t) is peaked near the edge x ~ w and it 
is of interest to estimate the width At of this peak in the 
velocity and in the voltage V{t) cx v{t)- The width At is 
definition dependent. For example, one can define it as the 
time interval between instants when v = Umax and time ri 
when V ~ (wmax + v)/2, where v is the background velocity 
due to the current /. In dimensionless units, v corresponds to 
dX/dr = 1. Thus we obtain 



6wXs + TT^ 
3wXs — TT^ 



with 



wXs 



< 1 



(47) 



(48) 



so that Ti < 0. Since |ri| ^ Xs ^ t^, we estimate the 
width of the velocity peak near the edge as At ^ ArXs/fi ^ 
Xg I /?, where the fraction of order unity in Eq. (|47i has been 
neglected. Therefore, the ratio of this width relative to the 
total crossing time tq of Eq. ( |43] | is 



At 

To 



Xs 



< 1. 



(49) 



IV. "COLD" AND "HOT" VORTEX CROSSINGS 

A vortex moving from the saddle point x — Xg the strip 
edge X ~ w, during the time tq = w/v, excites quasipar- 
ticles along its path by the mechanism described by Larkin 
and Ovchinnikovi^^ This mechanism is appropriate for dirty 
superconductors (for clean and intermediate clean regimes see 
Refs.l27landl28l). Since NbN films are inherently dirty, we can 
safely disregard the latter mechanism. We estimate the total 
energy transferred to quasiparticles during the time tq along 
the vortex path as 



e STT Ir 



N 



(50) 




This is, in fact, the work done by the Lorentz force on the 
vortex path of the length w — Xg- This energy is distributed 
near uniformly along the path at currents close to the critical 
current, because the vortex velocity varies weakly for most 
of the crossing, see Fig. 1 . In a belt of width £ along the y- 
axis with the volume Vb = £wd, the energy increase per unit 
volume is {87r^/e£){H^/8n){I/Ic). 

We now estimate the time tq. For a strip with resistivity 
p{Tc) = 240 /Ltficm, w = 120 nm, d = 4 nm, A = 45 ^m, 
and a bias current of the order Ic, the crossing time is roughly 
To ^ 10 ps and corresponding vortex speed is 12 km/s. This 
time is too short for any significant transfer of the electronic 
excitation energy into the substrate and surrounding strip area. 
Indeed, the phonon escape time was estimated as 160 ps in a 
strip of thickness d = 20 nm, whereas the electron-phonon 
relaxation time is about 17 psis^ During the time tq quasipar- 
ticles diffuse away from the vortex path by a short distance 
{DtoY^^ w 8 nm as estimated from the electronic specific 
heat Ce = 2.2 kJ/m'^K and the normal-state resistivity at 10 

Hence quasiparticles remain practically within the belt of 
volume Vb = £wd along the vortex path. The quasiparticle en- 
ergy density within the belt is {87r^/e£){H^ /87r){I/Ic). Tak- 
ing £ w 3^, we see that for I > I* ^ Ic/3 such an energy is 
sufficient to turn the belt normal causing a dark count in the 
photon detector. We call this process at high currents I > I* 
a "hot" vortex crossing. 

Therefore, we conclude that the superconducting strip with 
a bias current in the interval /* < J < Jc is unstable with 
respect to the transition into the normal state, that can be trig- 
gered by a vortex overcoming the barrier Clearly, photons 
can trigger such a transition as well. The photon efficiency in- 
creases as / approaches Ic and so does the rate of dark counts. 

In fact, the true critical current of a strip, below which the 
strip remains superconducting, is /*. At currents below /*, 
the superconducting state is stable, but remains resistive due 
to the presence of quasiparticles in normal cores of vortices 
crossing the strip. In this scenario a single vortex crossing 
leaves the strip in the superconducting state and thus we call 
this process a "cold" vortex crossing. 



FIG. 3. (Color online) The dark count rates of three SNSPDs at 5.5 K 
by Bartolf et al.— and fits based on Eqns. ( I27t and l l24t . The current 
is in units ly = ctic, where a = 0.72,0.77,0.60 for samples 1, 
2, and 3. Jc is the critical current defined as the current at which 
the barrier for vortex crossings vanishes. At low currents electronic 
noise in the measurement setup dominates over vortex crossings. 



V. COMPARISON WITH EXPERIMENTAL DATA 

In Fig. [3] experimental dark count rates are shown for three 
different NbN samples of SNSPDs,ii We fit the data using 
Eqs. ^} and by writing 



ln(i?/L) = ln(a 



ln[y($o//7ri^cr)], 

Q X 1/2 / 



(51) 
(52) 



The dimensions of samples 1, 2, and 3 are d = 6 nm, 
w = 53.4,82.9,170.6 nm, i = 73.9, 145.1, 141.4 ^m, re- 
spectively. The sheet resistance = 445,393,43157, and 
data were taken at T = 5.5 K. According to Bartolf et al., 
at low currents the data was dominated by electronic noise in 
the measurement circuit The data for samples 1 and 2 agree 
well with the theoretical results for high currents, while the 
data for sample 3 yield an unreasonably large exponent u. 

For sample 1 with fit parameter i' — $§/87r^A, we extract 
the Pearl length A(5.5K) — 57.1 fim, and from ln(a/L) we 
estimate the coherence length ^(5.5K) = 3.9 nm. The authors 
of Ref. [18| estimated ^(0) =4 nm from independent measure- 
ments of the upper critical field. They also estimated the Pearl 
length for NbN films of thickness d = 6 nm, A(0) = 65.1 
/im, from known resistivity p„ and the superconducting gap 
A(0) 2-3 meVi^ By using Eq. (O, we find the crit- 
ical current Ic = 20.1 pA defined as the current at which 
the energy barrier vanishes for vortex crossings. The authors 
of Ref. [la defined the "critical" current ly = 14.5 ^A us- 
ing the 1% voltage criterion (current at which resistance is 
1% of the normal one). We see that the critical current de- 
fined through such a voltage criterion is less than the critical 
current defined by the current at which the energy barrier van- 
ishes, ly ~ 0.72 Ic- 
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For the sample 2, we find the coherence length ^(5.5K) = 
4.33 nm and A(5.5K) = 51 /im. Independent estimates given 
in Ref.Hare ^(0) = 4.2 nm and A(0) = 59.2 firn; the critical 
current /c = 31.5 /i A, while /y = 0.77 /c^ '^ We conclude that 
our model for vortex crossing rates describes satisfactory the 
dark count rates in samples 1 and 2. 

Next we estimate the peak of the voltage pulse for / slightly 
below /*; 

Vpcak 7 , 

TTCAW 

and the duration of the pulse is Tpcak < 'I'o/cVpcak- For the 
sample 2 studied by Bartolf et alM we estimate V^poak ~ 0.8 
mV, while Tpoak ~ 3 ps slightly below /*. For comparison, 
dark counts are characterized by peak voltages of w 1 mV and 
by durations of several nanoseconds (FWHM ~ 2.5 nsi^). For 
dark counts the duration of pulses is caused by the current re- 
distribution and thus depends on the experimental setup used 
to detect the pulses. Note, that pulse duration differs signifi- 
cantly from that caused by single vortex crossing without for- 
mation of normal belt. 

The following experiment could, in principle, distinguish 
between regimes at / < /* and at / > /*: One induces 
a bias current in a thin-film ring and measures the magnetic 
flux in the ring as a function of time. For / > /*, a single 
vortex crossing destroys superconductivity and the flux van- 
ishes. The lifetime of this persistent current is 1/R and R is 
determined by Eg. (|27|i. If / < /*, the flux should decrease 
stepwise through multiple transitions between quantized cur- 
rent states /„, each transition corresponds to a single vortex 
crossing. In this case, the lifetime for the current /„ is l/i?„ 
where i?„ is given by Eq. ( |27l ) with / = /„. The total decay 
time of the initial current In will be r = ^n^- Foi" 

ID wires similar behavior due to phase slips was described by 
McCumber and Halperin^ 

In comparing theory and experiment, the issue of possible 
inhomogeneities of the thickness d and the width w is often 
raised. We note that the model developed here is only valid 
for w A. Each vortex in a narrow strip has mostly the 
kinetic energy of its supercurrents which are confined within 
an area of size ^ w x w. In other words, the model is not 
sensitive to inhomogeneities of d and of the edge roughness 
on scales small relative to w. 

Finally, it is worth mentioning that we assumed in this work 
that the strip temperature is equal to the bath temperature of 
the substrate. This may not always be the case in measure- 
ments of dark counts in photon detectors. After redistribu- 
tion of the bias current, the normal belt induced by a cross- 
ing vortex cools down. The strip can carry the supercon- 
ducting current equal to the bias current / only if the tem- 
perature drops below the value T* defined by the condition 
/* (T* ) = /. Slightly below T* vortices can cross the strip 
inside the warmer belt whose temperature is close to T* or 
inside the cooler areas whose temperature is that of the bath. 



The rate of vortex crossings is determined by both processes 
and the latter dominates only in the limit of very large L. 
Again, we emphasize that the measured rate is higher than the 
calculated rate, and the difference is larger for small currents 
because for them T* is higher 

VI. CONCLUSIONS 

In summary, we have found that the most plausible mecha- 
nism for dark counts in photon detectors is due to thermal fluc- 
tuations related to vortex crossings in the metastable current- 
carrying superconducting state, which is realized at bias cur- 
rents above some value /* ^ Ic/^- We conclude by listing 
our main results: 

(a) Vortices crossing the current-biased strip due to thermal 
fluctuations induce voltage pulses which can be detected ex- 
perimentally. The barrier for vortex crossings vanishes at the 
critical current defined by Eq. (IT4l . 

(b) In narrow and thin strips, the superconducting state is un- 
stable in the current interval I* < I < Ic and a transition into 
the normal state is triggered by vortices crossing the strip ac- 
companied by energy (heat) release. 

(c) We estimated the threshold for "hot" vortex crossings to 
be roughly /* w 7^/3. 

(d) Dark counts in current-biased superconducting strips re- 
ported in the literature were observed in the regime of 
metastable superconducting state. 

(e) At currents below /*, vortex crossings do not induce tran- 
sitions into the normal state, but still induce voltage pulses and 
the superconducting state is resistive due to the quasiparticles 
inside vortex cores of crossing vortices. We proposed a ring 
experiment, which allows to distinguish different decay pro- 
cesses of circular currents above and below /*. 

(f) We estimated the amplitude and duration of "cold" voltage 
pulses which can be detected below /*. 

Clearly it is desirable to test our theory by measuring I- 
V characteristics with a pulsed current technique to avoid 
heating. Further it will be interesting to study the rate and 
the shape of "cold" pulses at currents below /* at different 
temperatures and see their evolution from thermally induced 
crossings to quantum tunneling. 
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